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Abstract 

In an effort to examine the relationship between flare flux and corresponding 
CME mass, we temporally and spatially correlate all X-ray flares and CMEs in 
the LASCO and GOES archives from 1996 to 2006. We cross-reference 6 733 
CMEs having well-measured masses against 12 050 X-ray flares having position 
information as determined from their optical counterparts. For a given flare, we 
search in time for CMEs which occur 10-80 minutes afterward, and we further 
require the flare and CME to occur within ±45° in position angle on the solar 
disk. There are 826 CME/flare pairs which fit these criteria. Comparing the flare 
fluxes with CME masses of these paired events, we find CME mass increases with 
flare flux, following an approximately log-linear, broken relationship: in the limit 
of lower flare fluxes, log(CME mass) oc 0.68xlog(flare flux), and in the limit of 
higher flare fluxes, log(CME mass) oc 0.33xlog(flare flux). We show that this 
broken power-law, and in particular the flatter slope at higher flare fluxes, may 
be due to an observational bias against CMEs associated with the most energetic 
flares: halo CMEs. Correcting for this bias yields a single power-law relationship 
of the form log(CME mass) cx 0.70x log(flare flux). This function describes the 
relationship between CME mass and flare flux over at least 3 dex in flare flux, 
from «10~ 7 -10- 4 W mT 2 . 
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1. Introduction 



The physical connection between solar flares and coronal mass ejections has 
long been a topic of debate and ongoing research in solar physics. CMEs have 
been observed to occur in conjunction with flares and eruptive prominences 
(Mu nro et al, 1979{ |Webb and Hundhausen, 1987| ) and with helmet streamer 
disruptions ( Dryer, 1996 1. While these phenomena are not causally related, it 
can be definitively said that both X-ray flares and CMEs involve reconfiguration 
of com plex mag netic topologies within the corona via magnetic reconnection. 
Indeed. ISvest ka (2001) asserted that CMEs all share the same cause — magnetic 
field lines opening — and that the factor governing the properties of the resulting 
CME is the magnetic field strength in the region whence the CME originates. 
INindos and Andrewsl (|2004p proposed that helicity may be the link between 
flares and CMEs: CMEs have been observed to facilitate helicity loss, "carrying" 



high helicity magnetic flux from the Sun ( Chen et al, 1997| [Wood et al, 1999 



|Vourlidas et al, 2000[ ); |LaBonte, G corgoulis , and Ru st (2007) observed that ac- 
tive regions producing X-class flares generate enough helicity to match that lost 
via CME within hours to days after the flare. 

Efforts have been made to identify a physical link between the two phenom- 
ena by correlating properties of cospatial, contemporaneous flares and CMEs. 
Although flares themselves are not the cause of CMEs, their properties could 
serve as prediction tools for imminent CMEs and particle events, which are of 
much concern in the field of space weather. 

Statistical relationships between solar flares and CMEs are of interest to 
astronomers as large scale time-series X-ray observations have found solar-like 
X-ray activity on T Tauri Stars (hereafter referred to as TTS; solar-mass, pre- 
main sequence stars). Indeed, Haisch, Antunes, and Schmitt (1995) noted that 
stellar X-ray flare light curves behave similarly to solar X-ray flares, that is, 
the structure of the light curve (impulsive initial rise followed by exponential 
decay) is effectively identical, though the TTS flares are often several orders of 
magnitude more energetic. Based on these observations, theory and modeling 
of stellar flares has been guided by the premise that the physics behind both 
solar and stellar X-ray flares is the same (e.g., IReale et alX I1997[) . IPeres et al.\ 
(2001) demonstrated that stellar coronal X-ray emission could be reproduced 
by "scaling up" solar observations to include greater fractional coverage of the 
surface by active regions and solar-like coronal structures. Further underscoring 
the similarities in solar/stellar coronae, features similar to helmet streamers and 
slingshot prominences have indeed been found on TTS (e.o,. IMassi et aU 12008; 
ISkelly et al 

In many applications of TTS X-ray studies, the observable phenomenon is 
the X-ray flare, but the desired quantity is the associated mass loss. While 
solar CMEs do not shed large quantities of mass, young Suns («1 Myr) ex- 
hibit «3 orders of magnitude more energetic flares at a higher frequency, and 
the problem of how substantial angular momentum vis-a-vis mass is shed in 
these stars remains unresolved. Additionally, CME-like events on young stars 
could aid in understanding circumstellar disk evolution and planet formation: 
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for example, could CMEs on young stars be a mechanism for "Hash-heated" 
chondrule formation ( Miura and Nakamoto, 2 007)? 

To be clear, the motivation for this work is not to suggest a causal solar 
flare-CME relationship. We seek to quantify a relationship between solar flare 
fluxes and CME masses under the general premise that in some cases, flares 
and CMEs arise from common regions of complex magnetic topology and high 
field strength and thus may have correlated properties. Specifically, we aim to 
calibrate a relationship between solar flare flux and CME mass for the purpose 
of application to other stars possessing solar-like coronae. 



2. Archival Data 

We obtained archival measurements of flares and CMEs for the period of January 
1996 to December 2006. The LASCO CME database ( |Gopalswamy et ai, 2009l E| 
catalogs observations of the Large Angle and Spectrometric Coronagraph dating 
back to January, 1996. LASCO observations were complete for «83% of the 11 
years analyzed in this work, with 13 862 manually identified CMEs and their 
measured parameters from this period. Of particular interest in this work, 6 733 
CMEs with well-measured masses — i.e. not halo CMEs or CMEs which faded 
before entering the LASCO field of view — are reported (see Figure [T]), along 
with their linear speeds, accelerations, angular widths and central position angles 
(CPAs). CME masses are measured by converting brightness per pixel into mass 
per pixel; this calculation assumes the observed brightness is light Thomson 
scattered off of electrons in a mixture of 90% ionized Hydrogen and 10% Helium. 
All pixels within a manually defined CME shape are then integrated to obtain a 
total CME mass dVourlidas et ai, 200 0). The CME start time reported is when 
the CME is first detected in the C2 telescope field of view; C2 images the region 
from 2.0 to 6.0 R Q . 

Spanning more than three decades of observation, the Geostationary Opera- 
tional Environmental Satellite (GOES) flare databas^l reports the 1-8 A band 
full-disk X-ray flux at Earth; flare classifications are then applied based upon the 
peak flux in that bandpass. From 1996-2006, the database contains information 
for 22 674 flares. Of these, positions for optical counterparts are documented 
for approximately half of the flares in the database (Figure [5]). Of these 12 050 
flares, one A class, 3 638 B, 7 248 C, 1 056 M, and 107 X class flares are recorded. 



3. Determining Flare-CME Association 

As described above, we use the 6 733 CMEs from the LASCO CME catalog with 
well measured masses and the 12 050 flares with optical counterpart positions 
from the GOES flare database. 



1 This CME catalog is generated and maintained at the CDAW Data Center by NASA and The 
Catholic University of America in cooperation with the Naval Research Laboratory. SOHO is 
a project of international cooperation between ESA and NASA. 

z http://www.ngdc.noaa.gov/stp/SOLAR/ftpsolarflares.html 
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In determining whether a given flare and CME are associated, we utilize 
spatial and temporal criteria, requiring that both flare and CME occur within 
a set time window and angular separation. We initially set the time window 
generously to select CMEs which occur within ±2 h of a flare's start time (though 
we refine this below). Converting Stony hurst system flare positions (Cartesian) 
on the disk to spherical coordinates, we show in Figure[3]the angular separations 
of these time-matched CMEs and flares. A clear peak is seen about 0°± ~ 70°; 
we adopt a separation criterion of ±45° (e.g., both events occur in the same 
quadrant of the disk with respect to the CPA of the CME) as a value intermediate 
to the ±70° in our Figure [3] and the ks ±30° results of lYashiro et aL\ (|2"0Qg|l . 
Interestingly, but beyond the aims of this work, we note that there is a significant 
peak for flares which occur «180° in separation from the CME CPA, potentially 
indicative of large-scale correlated disruptions. 

After applying the ±2 h temporal and a ±45° angular separation cut, we 
assess whether the initially chosen, large time window can be narrowed (Figure 
|4]). Regardless of which flare time (start, peak or end) we use, there is a significant 
peak between 0-80 min in the resulting flare/CME pairs' time separations (CME 
start time - flare time). Other time separations surrounding these peaks appear to 
represent a "background" level at N «60 pairs per 10 min time bin; these could 
represent randomly matched flare-CME pairs which are not truly associated. 
Events appear to be the most strongly correlated when using the flare start time 
(solid, black histogram in Figure 3}; this is to say that the number of events 
in the 10-80 min bins are greater in number and there is less "background" 
variation than is seen for the other flare time choices. We therefore adopt the 
flare start time, but in any case this choice does not strongly affect our final 
results. 

It is important to note that the CME "start time" in fact corresponds not 
to the time when the CME was launched, but rather to the time the CME 
first appears in the C2 detector (at a height of «3.23 Rq on average). Using the 
linear speeds of the CMEs reported in the LASCO database, we can calculate the 
approximate CME travel times from the solar surface. The average travel time 
so computed is 75 ± 17 min, corresponding very well to the average time offset 
between the CME "start times" (detection times) and their associated flare times 
(Figure @| . Note that this simple calculation assumes negligible acceleration of 
the CMEs. Thus, while we cannot say with certainty that any one correlated 
flare/CME pair does not truly have a time delay, the data are broadly consistent 
with flare/CME start times that are simultaneous at the solar surface. 

Interpreting the 10-80 min peak of Figure 0] as a time offset region during 
which there is a higher probability of finding associated flares and CMEs, we 
narrow the time correlation window to accept CMEs which occur 10-80 min 
after a flare (vertical, triple dot-dashed lines, Figure SJ). This is consistent with 
the findings of lAndrews and Howard! (12001 jl . lMahrous et aZ.I (|2009)) ; in the latter 
work, the analysis used flare peak flux times and end times. This 10-80 min choice 
in temporal separation could potentially eliminate some true flare-CME pairs — 
previous studies (e.g., IHarrisonl 119911 Harrison! I1995[) have shown instances in 
which CMEs precede flares — but we do not consider these cases due to the 
absence of a distinct peak in the range of negative time offsets. The 10-80 min 
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window also evidently includes some "background" events. Potentially, these 
cases only represent a minority of flare-associated CMEs. 

Since we cannot tell from this data set whether a CME originates on the 
visible side of the limb or just beyond it, it is possible that some "back-sided" 
CMEs are spuriously associated with flares on the visible side. We can estimate 
the effect of back-sided CMEs using geometric arguments. First, we assume 
that CMEs observed by LASCO can be detected if they originate from a 260° 
longitudinal range about the Sun (i.e., from the full 180° of the visible disk, plus 
40° beyond both the eastern and western limbs), and we assume that CMEs are 
equally probable to originate from anywhere within this area. Thus, 33% of all 
CMEs (80°/240°) could be back-sided CMEs that are then potentially paired 
with flares on the visible solar disk. Our spatial cut requires ±45° coincidence 
in polar angle; assuming CMEs are equally likely to occur over all 180° of polar 
angle, this serves to eliminate half of all spurious correlations, reducing the 
number of back-sided CMEs potentially introduced to our correlated sample to 
16%. Further reducing this is our time-matching requirement, wherein we select 
only the CME-flare pairs occurring within 70 min, out of the full 240-min time 
window searched (i.e., 70/240 = 29%). Combining these, we potentially include 
5% of back-sided CMEs in our flare-CME pairs. Of course, some of these back- 
sided CMEs may in fact be truly associated with the observed flares; however, 
as a conservative statement, back-sided CMEs should represent at most a 5% 
contribution of spuriously flare-correlated CMEs in our sample and thus likely 
add only modestly to the scatter in our correlated quantities. 



4. Results 

4.1. Frequency of Associated Flares and CMEs 

With each of the above constraints applied, the total number of flare-CME pairs 
in the sample is winnowed as illustrated in Figure [5] The final sample contains 
826 flare-CME pairs, 737 of these are unique. As performed, the correlation 
allows for the selection of multiple CMEs per flare, and one CME may be 
found to be associated with multiple flares. With the data available, there is 
no way to distinguish which CME may "actually" be associated with the flare, 
or if indeed there should be only one associated. We therefore use all 826 pairs 
when comparing flare and CME properties, but we use only the 737 unique 
CMEs' properties when analyzing properties of flare-associated CMEs vs CMEs 
occurring without flares. 

Of the 6 733 CMEs we perform the correlation with, 737 unique CMEs are 
found to be associated with flares; put another way, we find 11% of CMEs 
to be flare-associated. Here we assess whether our applied criteria for corre- 
lation or observational effects could be decreasing the number of associations 
found, and in which variables these issues may create the most problems. In 
particular, the LASCO data gaps in addition to CME width and projection 
effects likely make our 11% an underestimate. For comparison, previous esti- 



mates of CMEs associated with flares are 40%, ( Munro et al, 1979 ) and 34% 
QSt. Cyr and Webb, 199l| . 
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4-.1.1. Data Gaps 

LASCO coverage was 83% complete from 1996-2006. We found that the following 
numbers of flares occurred during a LASCO data gap: 24 X, 119 M, 595 C, and 
803 B class. These flares could have been "lost" to our analysis, as the data gaps 
ended 10 min or more after the flare, within our 10-80 min temporal selection 
window. These values do not reflect, however, cases in which a CME did occur 
after the data gap ended and the flare was paired with that CME; we could not 
exclude these cases from the analysis as the technique does not allow for discern- 
ing which flares and CMEs are "actually" associated. Proportionately, more high 
energy flares occur during data gaps, likely because the particle events associated 
with the flares cause gaps in coverage by the LASCO C2 detector. B flares are 
excluded from this statement, as their energies are near the detection threshold 
for the GOES detector. An additional caveat to this is that the numbers of "lost" 
flares reported above do not account for simultaneity of flares; multiple flares 
could occur in a given data gap (thus, we cannot conclusively link a given class 
flare to a higher frequency of particle events), and we have counted each singly. 

4.1.2. Wide, Partial Halo, and Halo CMEs 

Surprisingly, we find that in the application of our time and angular separation 
constraints, most of the 107 X class flares are lost from the sample. Only seven 
of the 826 pairs have X class flares. This contradicts the results of l Andre ws 2003, 
who found an almost 100% X-flare-CME association^]. 

A limitation imposed by our data completeness requirement for the flares and 
CMEs analyzed in this work (see Section [5J we require flares to have positions, 
and CMEs to have well measured masses) is that we have excluded halo CMEs, 
as their masses cannot be measured via the same technique used for CMEs 
perpendicular to our line of sight. Additionally, our analysis excluded CMEs 
which had high errors in their mass measurements; these CMEs are illustrated 
in Figure |B] as a red histogram. Most of the CMEs with poorly constrained 
masses are what would be designated as "wide" or "partial halo" CMEs (i.e., 
with observed widths greater than 60°). 

In an effort to account for these selection effects in our analysis, we performed 
the correlation again, including wide, partial halo, and halo CMEs. This inclusion 
increased the CME sample size to 8 128 CMEs. Because they lack clearly defined 
CPAs, we did not require spatial coincidence for halo CMEs and flares, thereby 
including any halo CMEs which were time-correlated with flares. CMEs with 
poor mass measurements typically do have documented CPAs, so those events 
were required to be co-spatial with associated flares within ±45°, as previously. 
This resulted in 1 153 flare-CME pairs, with 53 X, 278 M, 608 C, and 214 B 
class flares. Of the 1 153 pairs, 1 028 of the CMEs are unique, resulting in a 13% 
flare-association; this is effectively indistinguishable from the association fraction 
observed when including only CMEs with good masses (Section 14. 1[) . While the 



3 Note, however, that Andrews 2003 set out specifically to define a sample of associated "big 
flares" and CMEs. 
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entire sample is increased by «40% and we don't find a greater fraction of CMEs 
to be associated with flares, the greatest percentage increases occur with higher 
flux flares. Including CMEs with poorly measured masses and halo CMEs results 
in «8 times as many X class flares, and almost double the M class flares. C and 
B class flares show only modest ~20% increases in representation in the sample. 
We assess below (Section 14.2.41) the main effect this could have on our results, 
and show a method to possibly correct for it. 

In summary, we lose «75% of X class flares to the combination of data gaps 
and observational selection effects, and this "loss fraction" is greatest for the 
strongest flares, decreasing with decreasing flare flux. Despite these factors po- 
tentially over-narrowing our correlated sample, we have demonstrated that the 
requisite time and position conditions applied here are sound in reasoning and 
represent a conservative approach which may only under-count flare-CME pairs. 

4.2. Correlations in the Matched Flare-CME Sample 

For the matched sample of 826 flare-CME pairs (see Section [3]), we assess po- 
tential relationships between five different flare and CME parameters: flare flux, 
CME linear speed, CME acceleration, CME width, and CME mass. Primarily, we 
seek to define a flare flux-CME mass relationship. We also perform additional 
correlations with our large flare-CME pair sample in order to compare with 
previous work and to elucidate the underlying physics. 

4-2.1. CME Linear Speed, Acceleration, and Flare Association 

Flare and eruptive prominence-associated CMEs have been observed to have 
higher linear speeds than those which are not associated with flares ( Gosling et ai, 1976| 



Moon et ai, 2003 ). To check whether our flare-CME pairs reflect this observed 
relationship, we show in Figure |7| the two distributions of CME linear speeds 
for our sample. The mean linear speed for CMEs associated with flares is 495±8 
km s ; CMEs not associated with flares have an average speed of 422±3 km 
s _1 (the error we report here is the standard deviation of the mean). A i-test 
( |Press et al, 1995] ) assesses the probability that two distributions have signif- 
icantly different means. The t statistic for these distributions is 8.5, and the 
corresponding null-hypothesis probability (that the two distributions have the 
same mean) is <CfCP 6 . 

Due to the asymmetric skew of the distributions, we verify this result by per- 
forming a non-parametric analog of the i-test, the Wilcoxon-Mann- Whitney test 
(also known as a U test). The statistic is Z, and the corresponding probability 
assesses the likelihood, when sampling two distributions, of preferentially finding 
a higher value in one distribution than another. The result of the U test corrob- 
orates the i-test: the Z statistic is -9.0, and the null-hypothesis probability is 
<ClO~ 6 . We conclude that these two distributions have statistically significantly 
different means, consistent with observations that flare-associated CMEs have 
higher linear speeds than CMEs not associated with flares. 

Figure [5] shows the relationship between CME speed and acceleration. Evi- 
dently, there is a trend of decreasing acceleration (or rather, increasing decel- 
eration since the mean accelerations are negative) with increasing CME linear 
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speed. To test this statistically, we perform a Spearman correlation test, finding 
the p coefficient to be —0.22 — a monotonically decreasing relationship — and the 
corresponding null-hypothesis probability <ClO -6 . In summary, we see that the 
fastest CMEs are associated with flares, and the fastest CMEs decelerate the 
most on average. 

4.2.2. CME Mass, Acceleration, Flare Flux 

In FigurelHlwe examine more closely the potential effect of CME mass on acceler- 
ation for CMEs associated with flares. The highest mass division includes CME 
masses >10 15 g, while the middle and lowest mass subsets span the ranges of 10 14 
g < CME mass < I0 15 g and CME mass < I0 14 g, respectively. Again utilizing 
the U test, we find that all three distributions have statistically indistinguishable 
means. 

In Figure [TOl CME acceleration in the flare-associated case and the converse 
scenario is shown. The mean of the CME-flare pair distribution is -1.8±0.1 m 
s~ 2 , while the mean of the distribution of CMEs without flares is 0.07±0.25 m 
s~ 2 (the errors on these quantities are the standard deviations of the mean). We 
use a i-test to determine whether the acceleration distributions have significantly 
different means. The t statistic is -2.5, the corresponding null-hypothesis proba- 
bility 0.01; this indicates these distributions do indeed have significantly different 
means. A U test was also performed; the Z statistic found to be 4.6, and the 
null-hypothesis probability is 1.9xl0 -6 . The CME-flare pairs have CMEs which 
tend to decelerate on average, while the CMEs not associated with flares are 
centered about acceleration. 

We next assess how the magnitude of this deceleration for flare-associated 
CMEs may relate to the flux of the associated flare. In Figurc[TTJ the relationship 
between CME acceleration and flare class is shown for the 826 CME-flare pairs. X 
and M flares are combined for this comparison, as there are only 7 X class flares in 
the 826 flare-CME paired sample. We employ the U test, making no assumption 
that these distributions are Gaussian. Comparing the mean accelerations of B 
class and weaker with C class flare-associated CMEs, the Z statistic is -1.9, 
with a corresponding null-hypothesis probability of 0.03. Applying the U test to 
intermediate strength flares, the difference in deceleration between C, M, and X 
class flare-associated CMEs is not statistically significant. However, testing the 
same group of B/A class flare-associated CMEs against M class and stronger 
flare-associated CMEs, the Z statistic is -2.8, its null-hypothesis probability 
0.002. Evidently, just as flare-associated CMEs decelerate on average, CMEs 
associated with the strongest flares show the most deceleration. 

4.2.3. CME Width, Mass, and Flare Flux 

We have so far shown above that higher energy flares are associated with faster 
CMEs that decelerate the most. To test whether higher energy flares are also 
associated with faster, more massive CMEs, we checked to see if the width of 
the CME was related to the flux of the associated flare. In Figure [TJJ we show 
that indeed, higher flux flares are associated with wider CMEs. With increasing 
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flare flux, the width of the associated CME increases: the mean widths (plus or 
minus the standard devations of the mean widths) of the flare-associated CMEs 
are 80±10°, 63±1.8°, 53±0.9°, and 42±1.4° for X, M, C, and B class flares, 
respectively. 

Considering projection effects, it is unsurprising then that as the mean CME 
width increases, so too does the occurrence of halo CMEs. Halo CME occurrence 
is also shown in Figure [T^l with 80% of X class flares being associated with halo 
CMEs, while M, C, and B class flares have halo CME associations of 21%, 6%, 
and 3%, respectively. A i-test shows the mean widths of CMEs associated with 
M, C, and B class flares are significantly different. Comparing the mean widths 
of CMEs associated with M and C class flares, the t statistic is 5.4, its null- 
hypothesis probability <ClO~ 6 . For the mean widths of CMEs associated with C 
and B class flares, the t statistic is 6.2, and the corresponding probability is also 
<1(T 6 . 

Finally, we find that as CME width increases, so too does the mass of the 
CME (Figure [T3"]) . Put another way, the widest CMEs aren't necessarily more 
diffuse than the more narrow CMEs. 

4.2.4. CME Mass, Flare Flux 

For the associated flares and CMEs, the relationship between CME mass and 
flare flux is shown in Figure [Ml Dividing the pairs by flare flux into four groups, 
we fit the resulting distribution with Gaussian functions. It is clear that the 
Gaussian centroids progress to higher CME mass as the flare flux sampled 
increases (Figure [T4]) . 

In Figure [151 the 826 flare/CME pairs are binned into fourteen equal sets 
of N — 59 pairs each. This division was chosen to maximize the number of 
events with which to perform statistics. Visually, it is apparent that the mean 
CME mass in each bin increases with flare flux. There is an apparent "knee" to 
the function around log(flare flux) « -5.5, so we use two linear regression fits 
to describe the functional relationship on both sides of this "knee." The first 
function is fit to bins 0-8, and the second to bins 7-13, counting from left to 
right. There is slight overlap in that precisely where the "knee" of the function 
occurs is unclear. The first linear fit is of the form: 

log(CME mass) = (18.5 ± 0.57) + (0.68 ± 0.10) x log(flare flux). (1) 

The second linear fit follows the form: 

log(CME mass) = (16.6 ± 1.30) + (0.33 ± 0.26) x log(flare flux). (2) 

We utilize Spearman's p correlation test to assess the degree of correlation of 
flare flux and CME mass. The test does not assume a particular form of corre- 
lation, rather, a p coefficient of unity indicates complete monotonic relationship 
of the variables tested. For the mean values in each bin, we find the p coefficient 
to be 0.99 with a null-hypothesis probability of <ClO -6 . This p coefficient nearly 
equal to one and the corresponding very low probability confirm that these are 
highly significantly correlated properties. 
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We also performed a Pearson linear correlation test. The Pearson coefficient 
ranges from —1 to +1 and evaluates how linearly related two parameter sets are. 
We obtain a value of r = 0.96 when comparing the mean masses and fluxes in 
Figure 1151 This strong correlation means one could nearly perfectly predict the 
mean behavior, but do note that there is a scatter about each mean mass of psI 
dex. 

To assess whether a broken fit is indeed necessary, we fit the correlated flux 
and mass bins of Figure[15]with a single linear fit. The xt value of this fit is 0.998, 
while the xt values for each of the fits above are 0.75 and 0.76, respectively. We 
also determine the robustness of the "knee" by re-binning the data and fitting 
again. The data divided this time into seven bins of N — 118, we fit a single line. 
This fit has a xt of 0.87. For comparison, we apply another broken log-linear fit 
to the data with the "knee" placed at the same point, log (flare flux)~ -5.5. Each 
of the functions better fits the data (xf, values of 0.16 and 0.07, respectively, 
potentially indicating we have "over- fit" the data). In both cases, we find that 
fitting the CME mass-flare flux relationship with a broken linear function indeed 
produces a better fit than simply fitting all of the data points with a single linear 
function. We also find that the "knee" is robust and does not depend upon how 
the data are binned. 

As mentioned in Section 14.1.21 observational effects could be causing "lost" 
flarc/CME pairs in the highest flare flux regime, which could in turn be causng 
the knee in the relationship of Figure [TS] In Section El. 1.21 we re- introduce wide, 
partial halo, and halo CMEs to the analysis. Halo CMEs lack mass measure- 
ments, but what we do know about halo CMEs allows us to assign an upper 
limit to possible masses. We have shown that the highest flux flares are most 
often associated with halo CMEs, and the highest flux flares tend to be associated 
with the widest CMEs (Figure [TJ]). We have also shown that the widest CMEs 
are the most massive (Figure ITS!) . Putting this expanded sample into a binned 
representation as in Figure [15l we obtain a range of CME masses for the given 
flare flux of the pair, and assign, as an upper bound, the maximum possible mass 
for that bin to the halo CME. In this way, we obtain an estimate of the greatest 
possible mean masses per bin when including halo CMEs. 

In Figure [T51 we reproduce Figure [TSl now including halo CMEs and CMEs 
with poorly constrained masses. The linear fit that describes this relationship is 
of the form: 

log(CME mass) = (18.67 ± 0.27) + (0.70 ± 0.05) x log(flare flux). (3) 

In this case, single power-law fit is sufficient to adequately describe the data. 
As this case assumes the maximum possible halo CME mass, the true CME 
mass/flare flux relationship likely lies in the parameter space bounded by the 
broken fit of Figure [T5] and the single fit shown in Figure [THl 

5. Discussion and Conclusions 

For flares and CMEs arising from the same active regions, one would antici- 
pate some of their properties to be correlated. For example, a less magnetically 
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complex active region may produce lower flux flares and lower mass CMEs. Con- 
versely, a highly magnetically complex active region may produce both higher 
energy flares and more massive CMEs. 

Our method for pairing flares and CMEs creates a very large data set with 
which to statistically examine whether flare and CME parameters are correlated. 
With this data set, we have found correlated properties of associated flares and 
CMEs which are consistent with previous observations, and we have furthermore 
shown a relationship between the flux of the flare and the mass of the CME. In 
summary, we have shown: 

1. CMEs associated with flares have higher average linear speeds than CMEs 
not associated with flares (Section 14.2.11 495±8 km s _1 and 422±3 km s _1 , 
respectively) , 

2. CMEs with higher linear speeds tend to decelerate; the relationship of accel- 
eration and linear speed is correlated, and monotonically decreasing (Section 

3. Flare-associated CMEs tend to decelerate: flare associated CMEs have a mean 
acceleration of -1.8±0.1 m s~ 2 while CMEs not associated with flares have a 
mean acceleration of 0.07±0.25 m s~ 2 (Section 14.2.21) . 

4. The magnitude of a flare-associated CME's deceleration increases with in- 
creasing flare flux: M class and stronger flare-associated CMEs have a mean 
acceleration of -3.7±1.3 m s~ 2 , while B class and weaker flare-associated 
CMEs have a mean acceleration of 0.14±0.8 m s -2 (Section [4.2.21) . 

5. A CME's mass is apparently unrelated to its acceleration (Section I4.2.2[) . 

6. The width of a flare-associated CME is directly correlated with the flux of 
the flare, with X class flare-associated CMEs being the widest (80°±10°) and 
B class flare-associated CMEs the most narrow (42°±1.4°), 

7. Wider CMEs are more massive than narrow CMEs (Section l4.2.31 Figure [T"3")) . 
and 

8. The stronger the associated flare is, the more massive the CME is (Section 
I4.2.4P : quantitatively, we find log(CME mass) oc 0.7x log(flare flux). 

The final bullet point in the list above could have very interesting implications 
for studies of young, solar-type stars, which exhibit high-frequency, high energy 
X-ray flaring activity. With the functional relationships illustrated in Figures 
[T5l and [TH1 it is now possible to extrapolate a CME mass/flare flux function 
into the range of flare fluxes observed for young stars (hatched region, Figure 
[T"5)) . Whether the second function of the broken log-linear relationship (Figure 
[T"5|) or the single function of the analysis including halo and poorly measured 
CME masses (Figure fl6|) would be used in this extrapolation remains to be seen; 
additional analyses of halo CME masses will be necessary to determine whether 
the "knee" exists for physical reasons and remains present when halo CMEs are 
included with better mass estimates. If the change in slope is physical, it could 
potentially represent a saturation point, approaching a set limit of possible CME 
mass. For young stars, it would be interesting to see whether this saturation point 
has any correlation with the saturation point in X-ray activity. 
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Figure 1. Summary of the LASCO CME database. From 1996 to 2006, there were 7 741 
CME mass measurements; 6 733 were well-constrained, and 1 008 were poorly constrained. 
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Figure 3. For flares and CMEs matched within a ±2 h time window, the angular separation 
of flares and CME central position angles shows a clear peak about 0° ±45° (vertical lines). 
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Figure 4. Assessment of which flare time to use for associating flares and CMEs; CT and FT 
refer to CME start time and flare time, respectively. We show here a ±2 h time constraint as 
well as the ±45° CME CPA-flare angle separation applied, the difference in each histogram 
being which flare time was used. Previous studies (e.g., Mahrous et al., 2009) found degree of 
flarc-CME association to be dependent upon whether flare start, peak, or end time was used 
in flare-CME pairing. Time offsets are shown here of CME start time minus flare start (solid, 
black), peak (dashed, red), and end times (dot-dashed, gray). We observe the strongest peak 
in number of pairs when using the flare start time. The peak of the black histogram is a factor 
of Ri3 above the "background" and well-defined in shape, spanning a range of time offsets from 
10-80 min (triple dot-dashed, light gray vertical lines). 



SOLA: cfcorr.tex; 2 November 2010; 0:22; p. 16 



Flare-CME Flux-Mass Correlation 




log(CME moss [g]) 

Figure 5. Change in the number and distribution of CME masses as flare/CME pairs are 
correlated. The white histogram bounded by dark gray shows the full initial distribution of 
well constrained CME mass measurements. Pairing flares with CMEs occurring 10-80 min after 
the flare start time, the data set is greatly reduced in number (filled, light gray histogram). The 
position criterion, CME and flare position angle equality within ±45°, leaves 826 flare/CME 
pairs (solid black line, white filled histogram). In red, we show CMEs not associated with flares 
by these criteria. 
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Figure 6. CME widths and frequency of "halo" CMEs. The distributions of widths for CMEs 
with well measured masses (solid, black line) and poorly constrained masses (dashed, red line) 
are shown. Masses are more difficult to measure, and thus the most poorly constrained, for the 
widest CMEs. All but three halo CMEs in our sample have no measured masses (dot-dashed, 
gray line). 
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Figure 7. As in Figure[5] the red and black distributions show CMEs not matched with flares 
and CME/flare pairs, respectively. We find that CMEs associated with flares have, on average, 
faster linear speeds (see Section |4.2,1| |, In this and subsequent figures, we report on the plot 
the mean of each distribution with the error on that mean (i.e. the standard deviation of the 
mean, cr/^/N). 
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Figure 8. For flare-associated CMEs, linear speeds and accelerations are shown as gray as- 
terisks. In bins 300 km s —1 wide in linear speed, mean speeds and accelerations are plotted as 
triangles (red). The figure has been cropped so the y-axis range of the triangles (red) can be 
clearly seen; 13 data points are not shown, lying beyond ±50 m s — 2 . While the magnitude of the 
slope of these points is rather small, its existence is significant: a Spearman's p correlation test 
shows these two quantities are highly significantly anticorrelated (p = -0.22, null-hypothesis 
probability -CIO" 6 ). We see here that the fastest CMEs tend to decelerate, while the slower 
CMEs accelerate through the C2 field of view (Section 14.2. Ill . 
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Figure 9. Relationship between CME mass and CME acceleration. The dot-dashed (red) 
histogram represents accelerations for the highest mass CMEs, with log(CME mass [g]) > 
15. The white histogram bounded by dark gray shows accelerations for CMEs with masses > 
10 14 g and <10 15 g. The lowest masses shown, <10 g, are represented by the solid (black) 
outlined histogram. There is no apparent dependence of CME acceleration on CME mass (see 
Section l4~Z2l) . 
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Figure 10. CME accelerations: CME/flare pairs (solid, black distribution) and CMEs not 
associated with flares (dot-dashed, red). The means of the two distributions (printed on the 
plot, with the associated standard deviations of each mean) are similar in magnitude but 
significantly different. The t statistic is -2.5, with a null-hypothesis probability of 0.01. A 
non-parametric U test confirms this, with a Z statistic of 4.6 and a null-hypothesis probability 
of 2xl0" 6 (see Section EH - 
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Figure 11. CME acceleration binned by flare flux. M and X class flares (log(flux [W m — 2 ]) > 
-5; combined here due to small number of X class flares in sample) are paired with CMEs that 
have accelerations as shown by the dot-dashed (red) distribution. C class flares (-6 < log(flux 
[W m -2 ]) < -5) have corresponding CME accelerations as shown by the white histogram 
bounded by dark gray. A and B class (log(flux [W m — 2 ]) < -6) flares' associated CMEs have 
accelerations as shown by the solid (black) histogram. We report the mean value of each 
distribution, with the standard deviation of each mean, on the figure. The B class flare/CME 
acceleration distribution has a mean significantly different than the higher energy flare asso- 
ciated CMEs; CMEs associated with lower energy flares appear to modestly accelerate, while 
those associated with higher energy flares decelerate (see Section l4.2.2ll . 
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Figure 12. Widths of CMEs associated with flares of a given flux. The mean widths for 
each flare class are printed on the plot with the standard deviations of each mean. The mean 
widths for CMEs associated with M, C, and B class flares are statistically significantly different 
from one another. Interestingly, 80% of CMEs associated with X class flares have "halo" 
designations, that percentage decreasing with decreasing flare flux (see Section l4.2.3ll . 
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Figure 13. CME width and mass. Each point represents a binned order of magnitude in mass, 
each abscissa value plotted at the bin's center. The diamonds (black) are the mean CME width 
for each mass bin, and triangles (blue) are median values. Error bars show the error on each 
mean value, while the gray shaded boxes show the standard deviation of each mean. As CME 
width increases, so does the mass of the CME (see Section |4. 2. 311 . 
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Figure 14. Grouping flare-CME pairs by flare flux, ascending flux from top to bottom, 
see a clear increase in the centroid value of CME mass (Section l4,2.4ll . 



SOLA: cfcorr.tex; 2 November 2010; 0:22; p. 26 



Flare-CME Flux-Mass Correlation 



CO 

o 



16. 



15.5 - 



1? 15.0 - 



14.5 - 



o 

^ 14.0 



15.5 - 



13.0 




-6 



-4 



-2 







log(Flore Flux [W m" 2 ]) 



Figure 15. Relationship between CME mass and flare flux. Black exes: the 826 flare/CME 
pairs are binned into equivalent boxes of N = 59. Each asterisk point is centered on the 
mean flux and mass per 59 pair bin. The flux error bars shows the minimum/maximum flare 
flux values spanned by that bin, and the mass error bars are the error on the mean (i.e., 
(T/y/W). Four red diamonds: abscissae are the mean flux values for each of the groups as set in 
Figure [141 ±x error is the standard deviation of that mean. The ordinates are the mass values 
corresponding to the peaks of the Gaussian fits in Figure 1141 their errors are the fit errors 
of the centroid. The light gray shaded boxes in the background are of arbitrary width, but 
show the standard deviation of the mean mass plotted in the foreground. For comparison, the 
hatched region represents the highest flux flares observed on young stars (see Section ^1 . Two 
linear functions are fit to the data, the first to bins 0-8, and the second to bins 7-13. There 
is slight overlap in that precisely where the "knee" of the function occurs is unclear. The first 
linear fit is of the form log(CME mass)=(18.5±0.57)+(0.68±0.10)xlog(flare flux). The second 
linear fit follows the form log(CME mass)=(16.6±1.30)+(0.33±0.26)xlog(flare flux). 
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Figure 16. As in Figure fT5l the relationship between flare flux and CME mass is shown. In 
this case, we include the sample of CMEs with poorly constrained mass measurements (Figure 
[l] dashed histogram), and see a relationship (diamonds, black) much like the one in Figure [T5l 
The triangles (red) indicate how the function would change were we to include halo CMEs; 
we assign halo CMEs the maximum mass from the bin occupied by its associated flare. This 
function, fit by the dot-dashed (black) line, represents the steepest slope we could anticipate 
assuming halo CMEs have masses which lie within the observed distribution. 
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